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Abstract 

We analyze charmonium and bottomonium production at fixed target experiments. We find 
that inclusion of color octet production channels removes large discrepancies between exper- 
iment and the predictions of the color singlet model for the total production cross section. 
Furthermore, including octet contributions accounts for the observed direct to total J /if) pro- 
duction ratio. As found earlier for photo-production of quarkonia, a fit to fixed target data 
requires smaller color octet matrix elements than those extracted from high-pt production at 
the Tevatron. We argue that this difference can be explained by systematic differences in the 
velocity expansion for collider and fixed-target predictions. While the color octet mechanism 
thus appears to be an essential part of a satisfactory description of fixed target data, impor- 
tant discrepancies remain for the Xci/Xc2 production ratio and J/ip (ip') polarization. These 
discrepancies, as well as, the differences between pion and proton induced collisions emphasize 
the need for including higher twist effects in addition to the color octet mechanism. 
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1 Introduction 



Quarkonium production has traditionally been calculated in the color singlet model (CSM) |L|. 
Although the model successfully describes the production rates for some quarkonium states, 
it has become clear that it fails to provide a theoretically and phenomenologically consistent 
picture of all production processes. In hadroproduction of charmonia at fixed target energies 
(y/s < 50 GeV), the ratio of the number of J/ip produced directly to those arising from decays 
of higher charmonium states is under-predicted by at least a factor five 0. The xa to Xc2 
production ratio is far too low, and the observation of essentially unpolarized J / if) and ip 1 can 
not be reproduced. At Tevatron collider energies, when fragmentation production dominates, 
the deficit of direct J /ip and ip' in the color singlet model is even larger. This deficit has been 
referred to as the '^'-anomaly' 0, Q. 

These discrepancies suggest that the color singlet model is too restrictive and that other 
production mechanisms are necessitated. Indeed, the CSM requires that the quark-antiquark 
pair that binds into a quarkonium state be produced on the time scale r ~ 1/tuq with the 
same color and angular momentum quantum numbers as the eventually formed quarkonium. 
Consequently, a hard gluon has to be emitted to produce a 3 Si state in the CSM and costs one 
power of a s /ir. Since the time scale for quarkonium formation is of order l/(mQV 2 ), where v 
is the relative quark-antiquark velocity in the quarkonium bound state, this suppression can 
be overcome if one allows for the possibility that the quark-antiquark pair is in any angular 
momentum or color state when produced on time scales r ~ 1/itlq. Subsequent evolution 
into the physical quarkonium state is mediated by emission of soft gluons with momenta of 
order rriQV 2 . Since the quark-antiquark pair is small in size, the emission of these gluons can 
be analyzed within a multipole expansion. A rigorous formulation [||] of this picture can be 
given in terms of non-relativistic QCD (NRQCD). Accordingly, the production cross section 
for a quarkonium state H in the process 

A + B — > H + X, (1) 

can be written as 

l 

(?h = J2j dx 1 dx 2 fi/A(x 1 )f j/B {x 2 ) -> H) , (2) 

i,j o 

^^)=E^ W ^' ( 3 ) 

n 

Here the first sum extends over all partons in the colliding hadrons and fa/ a etc. denote the 
corresponding distribution functions. The short- distance (x ~ S> l/(rriQv)) coefficients 

^QQ[n] describe the production of a quark-antiquark pair in a state n and have expansions 
in « s (2mg). The parameters^ (O^) describe the subsequent hadronization of the QQ pair 
into a jet containing the quarkonium H and light hadrons. These matrix elements can not be 
computed perturbatively, but their relative importance in powers of v can be estimated from 
the selection rules for multipole transitions. 

T Their precise definition is given in Sect. VI of §. 
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The color octet picture has led to the most plausible explanation of the ' ^'-anomaly' and 
the direct J/ip production deficit. In this picture gluons fragment into quark-antiquark pairs 
in a color-octet 3 Si state which then hadronizes into a ip' (or J/ip) @, 0, ||- Aside from 
this striking prediction, the color octet mechanism remains largely untested. Its verification 
now requires considering quarkonium production in other processes in order to demonstrate 
the process-independence (universality) of the production matrix elements (O^), which is an 
essential prediction of the factorization formula (|3|). 

Direct J/ ip and ip' production at large pt 3> 2rriQ (where tuq denotes the heavy quark mass) 
is rather unique in that a single term, proportional to {Og ( 3 S\)) , overwhelmingly dominates 
the sum (|3]). On the other hand, in quarkonium formation at moderate pt ~ 2mg at colliders 
and in photo-production or fixed target experiments (pt ~ lGeV), the signatures of color 
octet production are less dramatic, because they are not as enhanced by powers of n/a s or 
Pt/rriQ over the singlet mechanisms. Furthermore, theoretical predictions are parameterized 
by more unknown octet matrix elements and are afflicted by larger uncertainties. In particular, 
there are large uncertainties due to the increased sensitivity to the heavy quark mass close to 
threshold. (The production of a quark-antiquark pair close to threshold is favored by the rise 
of parton densities at small x.) These facts complicate establishing color octet mechanisms 
precisely in those processes where experimental data is most abundant. 

Cho and Leibovich |J studied direct quarkonium production at moderate pt at the Tevatron 
collider and were able to extract a value for a certain combination of unknown parameters 
(OgfSo)) and (£>f( 3 P )> (H = J/tp, i/j', T(1S), T(2S)). A first test of universality comes 
from photo-production [11], [EJ, where a different combination of these two matrix elements 
becomes important near the elastic peak at z ~ 1, where z = p ■ k^/p ■ k 7 , and p is the proton 
momentum. A fit to photo-production data requires much smaller matrix elements than those 
found in ||. Taken at face value, this comparison would imply failure of the universality 
assumption underlying the non-relativistic QCD approach. However, the extraction from 
photo-production should be regarded with caution since the NRQCD formalism describes 
inclusive quarkonium production only after sufficient smearing in z and is not applicable in 
the exclusive region close to z = 1, where diffractive quarkonium production is important. 

In this paper we investigate the universality of the color octet quarkonium production 
matrix elements in fixed target hadron collisions and re-evaluate the failures of the CSM 
in fixed target production after inclusion of color octet mechanisms. Some of the issues 
involved have already been addressed by Tang and Vanttinen |13[ and by Gupta and Sridhar 
fl~4| , but a complete survey is still missing. We also differ from [DJ in the treatment of polarized 
quarkonium production and the assessment of the importance of color octet contributions and 
from WM in the color octet short-distance coefficients. 



The paper is organized as follows: In Sect. 2 we compile the leading order color singlet 
and color octet contributions to the production rates for ip', xj, J /ip as well as bottomonium. 
In Sect. 3 we present our numerical results for proton and pion induced collisions. Sect. 4 is 
devoted to the treatment of polarized quarkonium production. As polarization remains one 



of the cleanest tests of octet quarkonium production at large pt Ju^, 16], we clarify in detail 
the conflicting treatments of polarized production in [|16|] and 0. Sect. 5 is dedicated to a 
comparison of the extracted color-octet matrix elements from fixed target experiments with 
those from photo-production and the Tevatron. We argue that kinematical effects and small- x 
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effects can bias the extraction of NRQCD matrix elements so that a fit to Tevatron data at 
large p t requires larger matrix elements than the fit to fixed target and photo-production data. 
The final section summarizes our conclusions. 



2 Quarkonium production cross sections at fixed target 
energies 

2.1 Cross sections 

We begin with the production cross section for ip' which does not receive contributions from 
radiative decays of higher charmonium states. The 2 — > 2 parton diagrams produce a quark- 
antiquark pair in a color-octet state or P-wave singlet state (not relevant to ip') and therefore 
contribute to ip' production at order a 2 v 1 . (For charmonium v 2 ~ 0.25 — 0.3, for bottomonium 
v 2 « 0.08 — 0.1.) The 2 — ► 3 parton processes contribute to the color singlet processes at order 
a 3 v 3 . Using the notation in (0): 



^)Tl 3 Oi 2 



(Of ('S„)) + ^j(0 8 *'( 3 P„)) + ^(Of( 3 P2» 



+ e(*i* 2 - 4m c 2 /,) (Ofi'SJ) z 2 



1 - z 2 + 2z\nz l-z 2 -2z\iaz 
+ 



(4) 



(1-z) 2 (1 + z) 3 

a(gq^ilj') = (5) 

= 27(2m!)^ 5(X1X2 " 4m " /g) (°* ,(3S ^ ^ 

Here z = (2m c ) 2 /(sxix 2 ), \/s is the center-of-mass energy and a s is normalized at the scale 
2m c . Corrections to these cross sections are suppressed by either a s /ir or v 2 . Note that the 
relativistic corrections to the color singlet cross section are substantial in specific kinematic 



regions z — > 0, 1 |L7|]. For y/s > 15 GeV these corrections affect the total cross section by 
less than 50% and decrease as the energy is raised Furthermore, notice that we have 
expressed the short-distance coefficients in terms of the charm quark mass, M^> ~ 2m c , rather 
than the true ip' mass. Although the difference is formally of higher order in v 2 , this choice is 
conceptually favored since the short- distance coefficients depend only on the physics prior to 
quarkonium formation. All quarkonium specific properties which can affect the cross section, 
such as quarkonium mass differences, are hidden in the matrix elements. 

The production of P-wave quarkonia differs from S"-waves since color singlet and color 
octet processes enter at the same order in v 2 as well as a s in general. An exception is xa, 
which can not be produced in 2 — > 2 parton reactions through gluon-gluon fusion in a color 
singlet state. Since at order a 2 , the xa would be produced only in a qq collision, we also 
include the gluon fusion diagrams at order a 3 , which are enhanced by the gluon distribution. 
We have for Xco, 

Hgg^xco) = ^^5(xiX2-4ml/s)±(orePo)) (7) 

3(2m Cj )' j s 
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for Xd, 



a(qq 



XcO) 

Xco) 



167T 3 « 2 



27(2m c ) 3 s 



5{x l x 2 -Am 1 Js) (O^SJ) 



(9) 



H99 -> Xcl) 



2ir 2 n 3 1 
" ~ e{x 1 x 2 -Am 2 J S )—{Or i ( 3 Pi)} 



X 



9(2m c ^ 5 

4z 2 In z (z 8 + 9z 7 + 26z 6 + 28z 5 + 17 z 4 + 7z 3 - 40z 2 - Az - 4 

(l + z) 5 (l-z) 4 

^ 9 + 39z 8 + 145z 7 + 251^ 6 + 119z 5 - 153z 4 - 17z 3 - U7z 2 - 8z + 10 



H9Q -> Xd) 
and for x c2 



3(1 -z) 3 (l + z) 



(10) 



81(2m c ) 5 

16vr 3 a 2 
27(2m c ) 3 s 



6(^x2 - 4m c 7s)^(Of 1 ( 3 P 1 )) 



mi 



—z In 2! + 



4^ 3 - 9z + 5 



5(x 1 x 2 -4m 2 / S )(0 8 ^ 1 ( 3 S 1 )) 



-» Xc2) 



7r 3 ry 2 . 

~ 5(x ia ; 2 -4m 2 / S ) — (0f 2 ( 3 P 2 )> 



1 



45(2m c ) 3 s 




(12) 
(13) 
(14) 



Note that in the NRQCD formalism the infrared sensitive contributions to the gg-induced 
color-singlet process at order a 3 are factorized into the color octet matrix elements (Og CJ ( 3 Si)) , 
so that the qq reactions at order a 3 are truly suppressed by a s . The production of P-wave 
states through octet quark-antiquark pairs in a state other than 3 Si is higher order in v 2 . 

Taking into account indirect production of J / ip from decays of ip' and Xcj states, the J / ip 
cross section is given by 



oj/v = cr(J/VW + Br (XcJ -» J/^X) a XcJ + Br(^' -> J/V^O ety' 

J=0,l,2 



(15) 



where 'Br' denotes the corresponding branching fraction and the direct J/ip production cross 
section <j(J/ip)dir differs from (see flD) only by the replacement of ip' matrix elements 
with J/ip matrix elements. Finally, we note that charmonium production through B decays 
is comparatively negligible at fixed target energies. 

The 2 — > 2 parton processes contribute only to quarkonium production at zero transverse 
momentum with respect to the beam axis. The transverse momentum distribution of H in 
reaction ([TJ) is not calculable in the p t < Aq C d region, but the total cross section (which 
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averages over all pt) is predicted even if the underlying parton process is strongly peaked at 
zero p t . 

The transcription of the above formulae to bottomonium production is straightforward. 
Since more bottomonium states exist below the open bottom threshold than for the charmo- 
nium system, a larger chain of cascade decays in the bottomonium system must be included. 
In particular, there is indirect evidence from T(3S') production both at the Tevatron [ITJ 



as 



well as in fixed target experiments (to be discussed below) that there exist yet unobserved 
X&(3P) states below threshold whose decay into lower bottomonium states should also be 
included. Our numerical results do not include indirect contributions from potential P-wave 
states below threshold. 

All color singlet cross sections compiled in this section have been taken from the review 
0. We have checked that the color octet short-distance coefficients agree with those given in 
||, but disagree with those that enter the numerical analysis of fixed target data in [|IJ[]. 

2.2 Matrix elements 

The number of independent matrix elements can be reduced by using the spin symmetry 
relations 



(Or j ( 3 Pj)) = (2J+l)(Or o ( 3 Po)) 

(OiCPj)) = (2J+l)(OtCP )) (16) 

(o 8 ^( 3 5 1 )) = (2j+i)(or( 3 5i)) 

and are accurate up to corrections of order v 2 (if) = J/ip,ip' - identical relations hold for 
bottomonium). This implies that at lowest order in a s , the matrix elements (O^^Sq)) and 
(Og( 3 Po)) enter fixed target production of J/ip and ip' only in the combination 

A 8 (H) = (O»CS )) + 4^( 3 P )> • (17) 

m Q 

Up to corrections in v 2 , all relevant color singlet production matrix elements are related to 
radial quarkonium wave functions at the origin and their derivatives by 

(Of ( s 50> = f |P(0)| 2 (Of ( 3 P )> = f |P'(0)| 2 . (18) 

We are then left with three non-perturbative parameters for the direct production of each 
S'-wave quarkonium and two parameters for P-states. 

The values for these parameters, which we will use below, are summarized in tables [TJ 
and |2|. Many of the octet matrix elements, especially for bottomonia, are not established 
and should be viewed as guesses. The numbers given in the tables are motivated as follows: 
All color singlet matrix elements are computed from the wavefunctions in the Buchmuller- 
Tye potential tabulated in [|lj| and using (|TJ|). Similar results within ±30% could be obtained 
from leptonic and hadronic decays of quarkonia for some of the states listed in the tables. The 
matrix elements (Og( 3 Si)) are taken from the fits to Tevatron data in |J with the exception 
of the 3S and 3P bottomonium states. In this case, we have chosen the numbers by (rather ad 
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ME 


J/v 




T(15) 


r(2S) 


T(35) 




1.16 
6.6 • 10" 3 
fitted 


0.76 
4.6 • lO- 3 
fitted 


9.28 
5.9- 10" 3 
5.0- 10" 2 


4.63 
4.1 • 10" 3 
3.0- 10~ 2 


3.54 
3.5- 10~ 3 
2.3 • 10~ 2 



Table 1: Matrix elements (ME) for the direct production of a S'-wave quarkonium H . All 
values in GeV 3 . 



ME 


XcO 


Xm(lP) 


Xbo(2P) 


Xw(3P) 


(0?( 3 Po))/m 2 Q 


4.4- 10~ 2 


8.5- 10~ 2 


9.9- 10- 2 


0.11 


W( 3 Si)> 


3.2 ■ 10~ 3 


0.42 


0.32 


0.25 



Table 2: Matrix elements (ME) for the direct production of a P-wave quarkonium H. All 
values in GeV 3 . 

hoc) extrapolation from the IS, 2S and IP, 2P states. The combination of matrix elements 
Ag(if) turns out to be the single most important parameter for direct production of J/ip and 
if)'. For this reason, we leave it as a parameter to be fitted and later compared with constraints 
available from Tevatron data. For bottomonia we adopt a different strategy since A 8 (H) is 
of no importance for the total (direct plus indirect) bottomonium cross section. We therefore 
fixed its value using the results of |J together with some assumption on the relative size of 
(Og(}S )) and ((9<f( 3 Po)) an d an ad hoc extrapolation for the 35* state. Setting A 8 (H) to 
zero for bottomonia would change the cross section by a negligible amount. 



3 Results 

Figs. | to | and table | summarize our results for the charmonium and bottomonium produc- 



tion cross sections. We use the CTEQ3 LO p(J parameterization for the parton distributions 



of the protons and the GRV LO f21| parameterization for pions. The quark masses are fixed 
to be m c = 1.5 GeV and m& = 4.9 GeV, as was done in |J. The strong coupling is evaluated 
at the scale /i = 2m,Q (Q = b, c) and chosen to coincide with the value implied by the param- 
eterization of the parton distributions (e.g., a s (2m c ) pa 0.23 for CTEQ3 LO). We comment 
on these parameter choices in the discussion below. The experimental data have been taken 
from the compilation in |T| with the addition of results from and the 800 GeV proton 
beam at Fermilab [fj3[ |24|. All data have been rescaled to the nuclear dependence A 0,92 for 
proton-nucleon collisions and A ' 87 for pion-nucleon collisions. All cross sections are given for 
xf > only (i.e. integrated over the forward direction in the cms frame where most of the 
data has been collected). 

3.1 if/ 

The total if)' production cross section in proton-nucleus collisions is shown in Fig. |I|. The 
color-singlet cross section is seen to be about a factor of two below the data and the fit, 
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10 20 30 40 



Figure 1: Total (solid) and singlet only (dotted) ip' production cross section in proton- nucleon 
collisions (xf > only). The solid line is obtained with Ag(ip') = 5.2 • 10~ 3 GeV 3 . 

including color octet processes, is obtained with 

A 8 (V>') = 5.2- l(T 3 GeV 3 . (19) 

The contribution from (Of ( 3 <Si)) is numerically irrelevant because gluon fusion dominates at 
all cms energies considered here. The relative magnitude of singlet and octet contributions is 
consistent with the naive scaling estimate 7r/a s -f 4 ~ 1 (The color singlet cross section acquires 
an additional suppression, because it vanishes close to threshold when 4m;?/ \x\X2s) — > 1). 

It is important to mention that the color singlet prediction has been expressed in terms 
of 2m c = 3 GeV and not the physical quarkonium mass. Choosing the quarkonium mass 
reduces the color singlet cross section by a factor of three compared to Fig. [I], leading to an 
apparent substantial ip' deficit|. As explained in Sect. 2, choosing quark masses is preferred 
but leads to large normalization uncertainties due to the poorly known charm quark mass, 
which could only be partially eliminated if the color singlet wave function were extracted from 
i/j' decays. If, as in open charm production, a small charm mass were preferred, the data could 
be reproduced even without a color octet contribution. Although this appears unlikely (see 
below), we conclude that the total ip' cross section alone does not provide convincing evidence 
for the color octet mechanism. If we neglect the color singlet contribution altogether, we obtain 
As(ip') < 1.0 • 10 -2 GeV 3 . This bound is strongly dependent on the value of m c . Varying m c 
between 1.3 GeV and 1.7 GeV changes the total cross section by roughly a factor of eight at 
t/s = 30 GeV and even more at smaller yfs. Compared to this normalization uncertainty, the 
variation with the choice of parton distribution and a s (fi) is negligible. This remark applies 
to all other charmonium cross sections considered in this section. 



■This together with a smaller value for the color singlet radial wavefvmction could at least partially explain 
the huge discrepancy between the CSM and the data that was reported in 




Figure 2: J/ip production cross sections in proton-nucleon collisions for xp > 0. The dotted 
line is the direct J/ijj production rate in the CSM and the dashed line includes the contribution 
from the color-octet processes. The total cross section (solid line) includes radiative feed-down 
from the Xcj and ip' states. The solid line is obtained with A 8 (J/ip) = 3.0 ■ 1CT 2 GeV 3 . 

3.2 J/ip 

The J/i/j production cross section in proton-nucleon collisions is displayed in Fig. A rea- 
sonable fit is obtained for 

A s (J/iP) = 3.0 ■ 10~ 2 GeV 3 . (20) 

We see that the color octet mechanism substantially enhances the direct J j ip production cross 
section compared to the CSM, as shown by the dashed and dotted lines in Fig. 0. The total 
cross section includes feed-down from Xcj states which is dominated by the color-singlet gluon 
fusion process. As expected from the cross section in Sect. 2, the largest indirect contribution 
originates from \c2 states, because Xci production is suppressed by one power of a s in the 
gluon fusion channel. The direct J/ip production fraction at y/s = 23.7 GeV (E = 300 GeV) is 



63%, in excellent agreement with the experimental value of 62% [55]. Note that this agreement 
is not a trivial consequence of fitting the color octet matrix element A 8 (J/ip) to reproduce 
the observed total cross section since the indirect contribution is dominated by color singlet 
mechanisms and the singlet matrix elements are fixed in terms of the wavefunctions of fl9|[ . 

One could ask whether the large sensitivity to the charm quark mass could be exploited to 
raise the direct production fraction in the CSM, thus obviating the need for octet contributions 
altogether? As shown in Fig. |3]this is not the case, since the charm mass dependence cancels in 
the direct-to-total production ratio. It should be mentioned, that expressing all cross sections 
in terms of the respective quarkonium masses increases (j(J/ip)dir/o'j/^, because M XcJ > Mjj^. 
However, the total color singlet cross section then decreases further and falls short of the data 
by about a factor five. We therefore consider the the combination of total J/ip production 
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Figure 3: Ratio of direct to total J/ip production in proton-nucleon collisions as a function of 
the charm quark mass in the CSM and after inclusion of color octet processes at E — 300 GeV. 
The experimental value is 0.62 ± 0.04. 



cross section and direct production ratio as convincing evidence for an essential role of color 
octet mechanisms for direct J/ip production also at fixed target energies. 

The comparison of theoretical predictions with the E705 experiment [^5| is summarized 
in Tab. [|. Including color octet production yields good agreement for direct J/\l/ production, 
as well as the relative contributions from all Xcj states and ip' '■ Note that the total cross 
section from p5[ is rather large in comparison with other data (see Fig. |2[). In the CSM, the 
direct production cross section of 7nb should be compared with the measured 89 nb, clearly 
demonstrating the presence of an additional numerically large production mechanism. Note 
also that our ip' cross section in the CSM is rather large in comparison with the direct J 





pN th. 


P N CSM 


piVexp. 


7i-iVth. 


7T"iVCSM 


7i iVexp. 




90 nb 


33 nb 


143 ± 21 nb 


98 nb 


38 nb 


178±21nb 


C(J r /ll))dir/vj/M> 


0.63 


0.21 


0.62 ±0.04 


0.64 


0.24 


0.56 ±0.03 


cr^>/cr{J/4>)dir 


0.25 


0.67 


0.21 ±0.05 


0.25 


0.66 


0.23 ±0.05 


X-fraction 


0.27 


0.69 


0.31 ±0.04 


0.28 


0.66 


0.37 ±0.03 


Xa/Xc2 ratio 


0.15 


0.08 




0.13 


0.11 


1.4 ±0.4 



Table 3: Comparison of quarkonium production cross sections in the color singlet model 
(CSM) and the NRQCD prediction (th.) with experiment at E = 300 GeV and E = 185 
GeV (last line only). The 'x-fraction' is defined by J2j=i,2^ T (Xcj — » J/ipX) cr XcJ /aj/^. The 
'Xci/x C 2 '-ratio is defined by Br(x c i -> J/ipX) a Xcl /(Bi(xc2 -> J/ipX) a Xc2 ). 
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Figure 4: Total (solid) and singlet only (dotted) ip' production cross section in pion-nucleon 
collisions (xp > only). The solid line is obtained with Ag(ip') = 5.2 • 10~ 3 GeV 3 . 



cross section in the CSM. A smaller value which compares more favorably with the data could 
be obtained if one expressed the cross section in terms of quarkonium masses ||. From the 
point of view presented here, this agreement appears coincidental since the cross sections are 
dominated by octet production. 

Perhaps the worst failure of the theory is the Xci to Xc2 ratio in the feed-down contribution 
that has been measured in the WAll experiment at E — 185 GeV |p6 |. We see that the 
prediction is far too small even after inclusion of color octet contributions. The low rate of 
Xi production is due to the fact, as already mentioned, that the gluon-gluon fusion channel 
is suppressed by a s /n compared to Xc2 due to angular momentum constraints. Together with 
J/i[) (and if/) polarization, discussed in Sect. 4, the failure to reproduce this ratio emphasizes 
the importance of yet other production mechanisms, presumably of higher twist, which are 
naively suppressed by Aqcd/^Ic 0j- 



3.3 Pion-induced collisions 

The ifi' and J/ip production cross section in pion-nucleon collisions are shown in Figs. |] and 
|5|. The discussion for proton-induced collisions applies with little modification to the pion 
case. A breakdown of contributions to the J / ip cross section at E — 300 GeV is given in 
table ^. The theoretical prediction is based on the values of A 8 (H) extracted from the proton 
data. Including color octet contributions can add little insight into the question of why the 
pion-induced cross sections appear to be systematically larger than expected. This issue has 
been extensively discussed in Q. The discrepancy may be an indication that, either the gluon 
distribution in the pion is not really understood (although using parameterizations different 
from GRV LO tends to yield rather lower theoretical predictions), or that a genuine difference 
in higher twist effects for the proton and the pion exists. 
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Figure 5: J/ip production cross sections in pion-nucleon collisions for xp > 0. Direct J/ip 
production in the CSM (dashed line) and after inclusion of color-octet processes (dotted line). 
The total cross section (solid line) includes radiative feed-down from the \ c j and ip' states. 
The solid line is obtained with A 8 (J/ip) = 3.0 • 1(T 2 GeV 3 . 

3.4 T(nS) 

If higher twist effects are important for fixed target charmonium production, their importance 
should decrease for bottomonium production and facilitate a test of color octet production. 
Unfortunately, data for bottomonium production at fixed target energies is sparse and does 
not allow us to complete this test. 

Due to the increase of the quark mass, bottomonium production differs in several ways 
from charmonium production, from a theoretical standpoint. The relative quark- ant iquark 
velocity squared decreases by a factor of three, thus, the color octet contributions to direct 
production of T(nS) are less important since they are suppressed by v 4 (at the same time 
a s {2mQ) decreases much less). The situation is exactly the opposite for the production of 
P-wave bottomonia. In this case the color singlet and octet contributions scale equally in 
v 2 . The increased quark mass, together with an increased relative importance of the octet 
matrix element (Og b0 ( 3 Si)) (extracted from Tevatron data in ||) as compared to the singlet 
wavefunction (compare Xco with Xbo in Tab. ||), leads to domination of quark-antiquark pair 
initiated processes. Consequently, the direct T(nS) production cross section is at least a factor 
ten below the indirect contributions from x^-decays. This observation leads to the conclusion 
that the number of T(3S) observed by the E772 experiment [E7| can only be explained if 



X&j(3P) states that have not yet been observed directly exist below the open bottom threshold. 
Such indirect evidence has also been obtained from bottomonium production at the Tevatron 
collider 



To obtain our numerical results shown in Fig. ^|, we assumed that these Xbj(3P) states 
decay into Y(3iS) with the same branching fractions as the corresponding n = 2 states. The 
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Figure 6: Total (direct plus indirect) T(nS) production cross sections (for xp > 0), consec- 
utively summed over n. The data point refers to the sum of n = 1, 2, 3. 



total cross sections are compared with the experimental value 195 ± 67pb/nucleon obtained 
from |23| at E = 800 GeV for the sum of T(nS), n = 1, 2, 3 and show very good agreement. 
The color-singlet processes alone would have led to a nine times smaller prediction at this 
energy. We should note, however, that integration of the a^-distribution for T(IS') production 
given in [27] indicates a cross section about two to three times smaller than the central value 



quoted by P^| . The theoretical prediction for the relative production rates of Y(IS') : T(2S') : 
T(3S) is 1 : 0.42 : 0.30 to be compared with the experimental ratio |27j 1 : 0.29 : 0.150. This 
comparison should not be over interpreted since it depends largely on the rather uncertain 
octet matrix elements for P-wave bottomonia. Due to lack of more data we also hesitate to 
use this comparison for a new determination of these matrix elements. 



4 and J/ip Polarization 

In this section, we deal with ip' and J/ip polarization at fixed target energies and at colliders 
at large transverse momentum. Before returning to fixed target production in Sect. 4.2, we 
digress on large-p t production. We recall that, at large pf ^> 4mg, ip' and direct J/ip produc- 
tion is dominated by gluon fragmentation into color octet quark-antiquark pairs and expected 
to yield transversely polarized quarkonia [[HJ. The reason for this is that a fragmenting gluon 
can be considered as on-shell and therefore transverse. Due to spin symmetry of NRQCD, 
the quarkonium inherits the transverse polarization up to corrections of order Am 2 c /p\ and v 4 . 
Furthermore, it has been shown |TIJ that including radiative corrections to gluon fragmenta- 



tion still leads to more than 90% transversely polarized ip' (direct J/tp). Thus, polarization 

§ These numbers were taken from the raw data with no concern regarding the differing efficiencies for the 
individual states. 
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provides one of the most significant tests for the color octet production mechanism at large 
transverse momenta. At moderate pf ~ 4mj?, non-fragmentation contributions proportional 
to (O^^So)) and {0$( 3 Pj)) are sizeable |J. Understanding their polarization yield quanti- 
tatively is very important since most of the p t -integrated data comes from the lower p r region. 
The calculation of the polarization yield has also been attempted in |J . However, the method 
used is at variance with jlBJ and leads to an incorrect result for S'-wave quarkonia produced 
through intermediate quark- ant iquark pairs in a color octet P-wave state. In the following 
subsection we expound on the method discussed in |TB[ and hope to clarify this difference. 



4.1 Polarized production 

For arguments sake, let us consider the production of a ip' m a polarization state A. This state 
can be reached through quark- ant iquark pairs in various spin and orbital angular momentum 
states, and we are led to consider the intermediate quark-antiquark pair as a coherent super- 
position of these states. Because of parity and charge conjugation symmetry, intermediate 
states with different spin S and angular momentum L can not interfere^, so that the only 
non-trivial situation occurs for 3 Pj-states, i.e. S — 1, L — 1. 

In it is assumed that intermediate states with different JJ Z , where J is total angular 
momentum do not interfere, so that the production cross section can be expressed as the sum 
over JJ Z of the amplitude squared for production of a color octet quark-antiquark pair in a 
3 Pjj z state times the amplitude squared for its transition into the ip'. The second factor can 
be inferred from spin symmetry to be a simple Clebsch- Gordon coefficient so that 

4 ] ~ E^WjJ) \(JWJm - A); 1A)| 2 . (21) 

JJz 

We will show that this equation is incompatible with spin symmetry which requires interference 
of intermediate states with different J. 

A simple check can be obtained by applying fl2"T| ) to the calculation of the gluon frag- 
mentation function into longitudinally polarized if}'. Since the fragmentation functions into 
quark-antiquark pairs in a 3 Pjj z state follow from by a change of color factor, the sum in 
(pl|) can be computed. The result not only differs from the fragmentation function obtained 
in |16| but contains an infrared divergence which can not be absorbed into another NRQCD 
matrix element. 



To see the failure of ( pl|) more clearly we return to the NRQCD factorization formalism. 
After Fierz rearrangement of color and spin indices as explained in ||, the cross section can 
be written as 



^ ~ H ai;hj ■ . (22) 

In this equation H ai -bj is the hard scattering cross section, and S a i-bj is the soft (non-perturba- 
tive) part that describes the 'hadronization' of the color octet quark pair into a ip' plus light 
hadrons. Note that the statement of factorization entailed in this equation occurs only on the 
cross section and not on the amplitude level. The indices ij and ab refer to spin and angular 

^Technically, this means that NRQCD matrix elements with an odd number of derivatives or spin matrices 
vanish if the quarkonium is a C or P eigenstate. 
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momentum in a Cartesian basis L a Si (a,i = 1,2,3 = x,y,z). Since spin-orbit coupling is 
suppressed by v 2 in the NRQCD Lagrangian, L z and S z are good quantum numbers. In the 
specific situation we are considering, the soft part is simply given by (the notation follows 

10) 

S% = (0\x^ t T A D a ) if) afa$ ^a,T A (-'- D b ) X |0> , (23) 

where Om destroys a if)' in an out-state with polarization A. To evaluate this matrix element at 
leading order in v 2 , we may use spin symmetry. Spin symmetry tells us that the spin of the if)' 
is aligned with the spin of the cc pair, so Sj£) b j oc e l *(A)e?(A). Now all vectors 5^.L can depend 
on have been utilized, and thus by rotational invariance, only the Kronecker symbol is left to 
tie up a and b. The overall normalization is determined by taking appropriate contractions, 
and we obtain 

S% = (OfCP ))5 oh e i *(X)e\X). (24) 

This decomposition tells us that to calculate the polarized production rate we should project 
the hard scattering amplitude onto states with definite S z = A and L z , square the amplitude, 
and then sum over L z {J2L z e a(L z )e (L z ) = <5 ab in the rest frame). In other words, the soft 
part is diagonal in the L Z S Z basis. 

It is straightforward to transform to the JJ Z basis. Since J z = L z + S z , there is no 
interference between intermediate states with different J z . To see this we write, in obvious 
notation, 

^ (A) ~ E , • S ?lj'Ji . (25) 

and using (p4|) obtain, 

S?Ij,j> = {Oi\'P Q ))Y.^M ] \X\JJ z ){J'J' z \lM-l\), (26) 

M 

which is diagonal in (JJ Z )(J'J' Z ) only after summation over A (unpolarized production). In 
general, the off-diagonal matrix elements cause interference of the following JJ Z states: 00 
with 20, 11 with 21 and 1(— 1) with 2(— 1). While the diagonal elements agree with (pl[), the 
off-diagonal ones are missed in (|2T|). 

To assess the degree of transverse if)' (direct J/ip) polarization at moderate p t , the calcu- 
lation of should be redone with the correct angular momentum projections. 



4.2 Polarization in fixed target experiments 

Polarization measurements have been performed for both if) [^| and if)' p9f production in 
pion scattering fixed target experiments. Both experiments observe an essentially flat angular 
distribution in the decay if) — > (if) = J/ip, if)'), 

, da - <xl + a cos 2 9 , (27) 
a cos 
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where the angle 6 is defined as the angle between the three-momentum vector of the positively 
charged muon and the beam axis in the rest frame of the quarkonium. The observed values for 
a are 0.02 ±0.14 for if)', measured at y/s = 21.8 GeV in the region x F > 0.25 and 0.028 ±0.004 
for J /if) measured at y/s = 15.3 GeV in the region xf > 0. In the CSM, the J/ip's are predicted 
to be significantly transversely polarized 0, in conflict with experiment. 

The polarization yield of color octet processes can be calculated along the lines of the 
previous subsection. We first concentrate on if)' production and define £ as the fraction of 
longitudinally polarized if)' . It is related to a by 

« = W" (28) 

For the different intermediate quark-antiquark states we find the following ratios of longitu- 
dinal to transverse quarkonia: 



(29) 



where the number for the singlet process (first line) has been taken from ||[§. Let us a dd the 
following remarks: 

(i) The ^j^-subprocess yields pure transverse polarization. Its contribution to the total 
polarization is not large, because gluon-gluon fusion dominates the total rate. 

(ii) For the 3 Pj 8 ' ) -subprocess J is not specified, because interference between intermediate 
states with different J could occur as discussed in the previous subsection. As it turns out, 
interference does in fact not occur at leading order in a s , because the only non- vanishing 
short- distance amplitudes in the JJ Z basis are 00, 22 and 2(— 2), which do not interfere. 

(iii) The 1 S'o 8 ^-subprocess yields unpolarized quarkonia. This follows from the fact that 
the NRQCD matrix element is 

(0tfT A afa$ ^T A X \0) = i (Of (%)> , (30) 



3 sf 


1 : 3.35 


i = 


0.23 


lo(8) 
D 


1 : 2 


£ = 


1/3 


3p(8) 


1 : 6 


t = 


1/7 


3g(8) 


: 1 




= 



independent of the helicity state A. At this point, we differ from [13], who assume that this 



channel results in pure transverse polarization, because the gluon in the chromomagnetic 
dipole transition 1 S , q 8 ' ) — > ^Sf^ + g is assumed to be transverse. However, one should keep in 
mind that the soft gluon is off-shell and interacts with other partons with unit probability prior 
to hadronization. The NRQCD formalism applies only to inclusive quarkonium production. 
Eq. ( |3~0D then follows from rotational invariance. 

(iv) Since the 3 -Pj 8 ' ) and ^o^-subprocesses give different longitudinal polarization fractions, 
the if)' polarization depends on a combination of the matrix elements (Of ( 1 5'o)) and (Of ( 3 Po)) 
which is different from A s (ip'). 

To obtain the total polarization the various subprocesses have to be weighted by their 
partial cross sections. We define 

"This number is ^F-dependent and we have approximated it by a constant at low Xp, where the bulk data 
is obtained from. The polarization fractions for the octet 2 — ► 2 parton processes are x^-independent. 



15 



Ss(H) 



and obtain 



(31) 



0.23 



+ 



0.16 + 0.11 5 8 (ip'), 



ay 



(32) 



where the last line holds at y/s = 21.8 GeV (The energy dependence is mild and the above 
formula can be used with little error even at y/s = 40 GeV). Since < S 8 (H) < 1, we have 
0.16 < f < 0.27 and therefore 



0.15<a<0.44. (33) 

In quoting this range we do not attempt an estimate of 5g(if)'). Note that taking the Teva- 
tron and fixed target extractions of certain (and different) combinations of (0 8 ( 1 5'o)) and 
(Og ( 3 -Po)) seriously (see Sect. 5.1), a large value of S 8 (ip') and therefore low a would be 
favored. Within large errors, such a scenario could be considered consistent with the mea- 
surement quoted earlier. From a theoretical point of view, however, the numerical violation 
of velocity counting rules implied by this scenario would be rather disturbing. 

In contrast, the more accurate measurement of polarization for J /if) leads to a clear dis- 
crepancy with theory. In this case, we have to incorporate the polarization inherited from 
decays of the higher charmonium states Xcj and if)'. This task is simplified by observing that 
the contribution from Xco and xa feed-down is (theoretically) small as is the octet contribu- 
tion to the Xc2 production cross section. On the other hand, the gluon-gluon fusion process 
produces Xci states only in a helicity ±2 level, so that the J /if) in the subsequent radiative 
decay is completely transversely polarized. Weighting all subprocesses by their partial cross 
section and neglecting the small if}' feed-down, we arrive at 



£ = 0.10 + 0.11 5 8 (J/if>) (34) 

at y/s = 15.3 GeV, again with mild energy dependence. This translates into sizeable transverse 
polarization 



0.31<a< 0.63. (35) 

The discrepancy with data could be ameliorated if the observed number of Xci from feed-down 
were used instead of the theoretical value. However, we do not know the polarization yield of 
whatever mechanism is responsible for copious Xci production. 

Thus, color octet mechanisms do not help to solve the polarization problem and one has to 
invoke a significant higher-twist contribution as discussed in 0j . To our knowledge, no specific 
mechanism has yet been proposed that would yield predominantly longitudinally polarized if)' 
and J /if> in the low xp region which dominates the total production cross section. One might 
speculate that both the low Xcv/Xci ratio and the large transverse polarization follow from 
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the assumption of transverse gluons in the gluon-gluon fusion process, as inherent to the 
leading-twist approximation. If gluons in the proton and pion have large intrinsic transverse 
momentum, as suggested by the prspectrum in open charm production, one would be naturally 
led to higher-twist effects that obviate the helicity constraint on on-shell gluons. 



5 Other processes 

Direct J/ip and i/j' production is sensitive to the color octet matrix element A 8 (H) defined in 
(|H|). In this section we compare our extraction of A 8 (if) with constraints from quarkonium 
production at the Tevatron and in photo-production at fixed target experiments and HERA. 

5.1 Quarkonium production at large p t 

An extensive analysis of charmonium production data at pt > 5 GeV has been carried out 
by Cho and Leibovich || M, who relaxed the fragmentation approximation employed earlier 
H 0| . At the lower p t boundary, the theoretical prediction is dominated by the and 3 Pj 8 ^ 
subprocesses and the fit yields 



<(9 8 J/ VSo)> + —.(O^CPo)) = 6.6 • 10- 2 
(Oi'CS )) + ^{Oi'( 3 P )) =1.8-10~ 2 , (36) 



to be compared with the fixed target valuesQ 



(O^fSo)) + ^(O J 8 /lp ( 3 P )) = 3.0 • 10- 2 
Tnr c 

(Ot'CS )) + — 2 (Ol'CP )) =0.5-10- 2 . (37) 

If we assume (e> 8 J/ YSo)> = (O J 8 ^( 3 P )}/m 2 c , the fixed target values are a factor seven (four) 
smaller than the Tevatron values for J ftp (ip'). The discrepancy would be lower for the radical 

choice (e> 8 J/V, ( 3 P )> = 0. 

While this comparison looks like a flagrant violation of the supposed process-independence 
of NRQCD production matrix elements, there are at least two possibilities that could lead to 
systematic differences: 

(i) The 2 — > 2 color octet parton processes are schematically of the form 

m.±, S(xiX2 s-M}), (38) 



** Since there is a strong correlation between the charm quark mass and the extracted NRQCD matrix 
elements, we emphasize that both (|6|) and ( p^ ) as well as (^) below have been obtained with the same 
m c = 1.5 GeV (or m c — 1.48 GeV, to be precise). On the other hand, the apparent agreement of predictions 
for fixed target experiments with data claimed in |14| is obtained from (p9|) in conjunction with m c = 1.7 GeV. 
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where Mf denotes the final state invariant mass. To leading order in v 2 , we have Mf = 2m c . 
Note, however, that this is physically unrealistic. Since color must be emitted from the quark 
pair in the octet state and neutralized by final-state interactions, the final state is a quarkonium 
accompanied by light hadrons with invariant mass squared of order M 2 m (Mh + M^v 2 ) 2 
since the soft gluon emission carries an energy of order M^v 2 , where M# is the quarkonium 
mass. The kinematic effect of this difference in invariant mass is very large since the gluon 
distribution rises steeply at small x and reduces the cross section by at least a factor two. 
The 'true' matrix elements would therefore be larger than those extracted from fixed target 
experiments at leading order in NRQCD. Since the if)' is heavier than the J /if), the effect is more 
pronounced for if)' , consistent with the larger disagreement with the Tevatron extraction for if)'. 
Note that the effect is absent for large-pt production, since in this case, X1X2S > 4pf 3> Mj. If 
we write Mf = 2m c + 0(v 2 ), then the difference between fixed target and large-pt production 
stems from different behaviors of the velocity expansion in the two cases. 

(ii) It is known that small-x effects increase the open bottom production cross section 
at the Tevatron as compared to collisions at lower y/s. Since even at large p t , the typical 
x is smaller at the Tevatron than in fixed target experiments, this effect would enhance the 
Tevatron prediction more than the fixed target prediction. The 'true' matrix elements would 
therefore be smaller than those extracted from the Tevatron in ||. 

While a combination of both effects could well account for the apparently different NRQCD 
matrix elements, one must keep in mind that we have reason to suspect important higher 
twist effects for charmonium production at fixed target energies. Theoretical predictions for 
fixed target production are intrinsically less accurate than at large p t , where higher-twist 
contributions due to the initial hadrons are expected to be suppressed by Ag C D/Pt (if not 
^QCD/Pt) rather than A QCD /m c . 



5.2 Photo-production 

A comparison of photo-production with fixed target production is more direct since the same 
combination of NRQCD matrix elements is probed and the kinematics is similar. All analyses 
lllOl, ITTl O find a substantial overestimate of the cross section if the octet matrix elements of 



(pq) are used. The authors of [[Tl[] fit 

(Oi /lp CS )) + -^(O^CPo)) = 2.0 • 10- 2 , (39) 
consistent with ( |3TD within errors, which we have not specified. While this agreement is 



reassuring, it might also be partly accidental since the extraction of is performed on the 
elastic peak, which is not described by NRQCD. Color octet mechanisms do not leave a clear 
signature in the total inelastic photo-production cross section. The authors of [[UJ argue that 
the color-octet contributions to the energy spectrum of J /if) are in conflict with the observed 
energy dependence in the endpoint region z > 0.7, where z = Ej/^/E^ in the proton rest 
frame. This discrepancy would largely disappear if the smaller matrix element of (37) or 



were used rather than (|36|). Furthermore, since in a color octet process soft gluons with energy 
Mhv 2 must be emitted, but are kinematically not accounted for, the NRQCD-prediction for 
the energy distribution should be smeared over an interval of size 8z ~ v 2 ~ 0.3, making the 
steep rise of the energy distribution close to z = 1 is not necessarily physical. 
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6 Conclusion 



We have reanalyzed charmonium production data from fixed target experiments, including 
color octet production mechanisms. Our conclusion is twofold: On one hand, the inclusion of 
color octet processes allows us to reproduce the overall normalization of the total production 
cross section with color octet matrix elements of the expected size (if not somewhat smaller) 
without having to invoke small values of the charm quark mass. This was found to be true 
for bottomonium as well as for charmonium. Comparing the theoretical predictions within 
this framework with the data implies the existence of additional bottomonium states below 
threshold which have not yet been seen directly. 

On the other hand, the present picture of charmonium production at fixed target energies 
is far from perfect. The Xci/Xc2 production ratio remains almost an order of magnitude too 
low, and the transverse polarization fraction of the J/ip and ip' is too large. We thus confirm 
the expectation of |2| that higher twist effects must be substantial even after including the 
octet mechanism. 

The uncertainties in the theoretical prediction at fixed target energies are substantial and 
preclude a straightforward test of universality of color octet matrix elements by comparison 
with quarkonium production at large transverse momentum. We have argued that small-x, 
as well as kinematic effects, could bias the extraction of these matrix elements in different 
directions at fixed target and collider energies. The large uncertainties involved, especially 
due to the charm quark mass, could hardly be eliminated by a laborious calculation of a s - 
corrections to the production processes considered here. To more firmly establish existence of 
the octet mechanism there are several experimental measurements which need to be performed. 
Data on polarization is presently only available for charmonium production in pion-induced 
collisions. A measurement of polarization at large transverse momentum or for bottomonium 
is of crucial importance, because higher twist effects should be suppressed. Furthermore, a 
measurement of direct and indirect production fractions in the bottom system would provide 
further confirmation of the color octet picture and constrain the color octet matrix elements 
for bottomonium. 
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